Kondo "underscreening" cloud: spin-spin correlations around a partially screened 

magnetic impurity 
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We consider the spatial spin correlations around a partially screened spin-1 magnetic moment 
in a metal exhibiting the underscreened Kondo effect. We find that the underscreening of the 
impurity spin results in spatial spin correlations that are more pronounced as compared to the 
fully screened Kondo effect; their power-law decay is weaker because of characteristic logarithmic 
corrections at large distances. The spin correlator also changes sign as a function of distance to the 
impurity allowing for ferromagnetic correlations between conduction electron spin density and the 
local moment. The numerical findings are shown to be in agreement with the predictions deriving 
from an effective ferromagnetic Kondo Hamiltonian. 
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Introduction — Spatial spin correlations around Kondo 
impurities - the Kondo spin screening cloud - and its pos- 
sible observability in experimental setups have been stud- 
ied intensively over the last decades [J, 0, [E 0, IE IE • 
Most of the work focused on the properties of spin corre- 
lations within the spin-i Kondo model, which is thought 
to capture the essential low-energy physics of magnetic 
ions in metals. It exhibits complete spin screening and 
Fermi liquid behavior below the Kondo temperature Tk, 
two characteristic features that also determine the long- 
distance behavior of the spin correlations. However, most 
transition metal and rare earth atoms have a spin mag- 
netic moment S greater than i in the ground state. 
Hence, when immersed in a metal they can potentially 
form underscreened Kondo impurity systems, even if the 
local moment degeneracy is partially lifted due to crystal 
field or spin-orbit splitting. 

The ground state of the spin-1 single-channel Kondo 
model - the prototype of underscreened Kondo models 
- is known to be effectively Nozieres Fermi liquid with a 
decoupled impurity spin-i. Indeed, a careful analysis of 
the conduction electron scattering matrix [E 0] revealed 
that the spectrum exhibits a one-to-one correspondence 
to that of a free Fermi gas. However, the residual ferro- 
magnetic interaction of the Nozieres quasiparticles with 
the residual impurity spin-^ results in non-analytic Fermi 
liquid corrections. Although there are well-defined quasi- 
particles at the Fermi surface, their decay rate vanishes 
with energy in a singular fashion. Correspondingly, the 
low-energy properties of the underscreened Kondo mod- 
els were termed as singular Fermi liquid behavior 
in particular, to contrast it with the genuine non-Fermi 
liquid characteristics of the overscreened Kondo mod- 
els. The anomalous corrections should affect all physical 
quantities that involve the conduction electron scattering 
matrix like, e.g., conductivity or, as we show below, the 
long-distance tail of the spin correlations. 



Generally, underscreened Kondo physics emerges if the 
effective number of screening channels that couple to the 
local moment carrying a spin S is smaller than 2S. Al- 
though the ground state of a Kondo spin in experimental 
systems is a fully screened singlet in most of the cases, 
there might be an extended temperature regime with an 
underscreened moment as the coupling of some of the 
channels can be very weak. For example, the screening of 
a magnetic impurity atom with spin S > \ adsorbed on 
metallic surfaces is restricted due to the lower symmetry 
near the surface potentially giving rise to underscreening 
ToL 11 1 . A realization in quantum dot setups was also 
proposed IE HI- If a quantum dot with a spin 5 = 1 
ground state is asymmetrically attached to leads, the spin 
is quenched in a two-stage process. The stronger cou- 
pled channel first screens half of the spin-1 at a temper- 
ature Tki , and the weaker coupled channel subsequently 
screens the residual spin-i at a lower temperature Tk2- 
For intermediate temperatures, T^ 2 < T < Tki, the 
quantum dot system is governed by the characteristics 
of the underscreened spin-1 Kondo model, and its singu- 
lar dynamics might be experimentally accessible in this 
temperature range. 

In this work, we focus on the spatial, equal-time spin- 
correlations of the conduction electrons around a spin-1 
impurity, that provide a snapshot of the screening cloud. 
In general, the characteristics of the screening cloud in 
Kondo models are less well studied than thermodynamic 
or transport properties. In the fully screened spin-i 
model, a number of theoretical results^, [E IE] support 
the expectation that the spatial spin-correlations exhibit 
a crossover at a length scale £x = Hvf/^bTk), where Vf 
is the Fermi velocity and Tk is the Kondo temperature [3|; 
for a review see Section 9.6 of Ref. [l[. Generally at zero 
temperature, spin-density correlations decay as a power- 
law with increasing distance x from the impurity. For the 
fully screened spin-^ impurity, the crossover is accompa- 
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nied by a change in the power-law exponent; while at 
short distances, x < £k , spin correlations decay as 1/x , 
with d the dimensionality of the conduction electron sys- 
tem, there is a stronger l/x d+1 decay at large distances, 
x > £,k 0, 0, 0] • Recently, it was shown that even the 
charge-density oscillations around the fully screened im- 
purity bear signatures of this crossover [7j, albeit with 
much weaker characteristics. On the experimental side, 
the detection of this crossover in spatial correlations has 
proved to be elusive so far. The most promising sys- 
tem for its observation is probably a magnetic impurity 
placed on a metallic surface subject to an STM study [7|. 

Comparatively, only little is known about the spatial 
correlations around a partially screened magnetic impu- 
rity. The question arises as to how the singular Fermi liq- 
uid properties of the underscreened spin-1 Kondo model 
modify the crossover at The scattering of the elec- 
trons off the unscreened residual moment will influence 
the long-distance behavior of spatial correlations, and the 
conduction electron cloud around a partially screened im- 
purity — the "underscreening" cloud — will thus differ 
from its fully screened counterpart. We study this ques- 
tion numerically with the help of an extention of Wil- 
son's numerical renormalization group (NRG) technique 
Q . We find that within the numerical accuracy the enve- 
lope of the equal-time impurity-spin - electron-spin cor- 
relator is a universal function of x/£k- In contrast to 
the spin-i model, its power-law decay obtains now only 
logarithmic corrections for x > £k> an d its long-distance 
tail is thus more pronounced than in the fully screened 
case. We demonstrate that this numerically observed tail 
of the "underscreening cloud" is in agreement with the 
predictions of an effective ferromagnetic spin-i Kondo 
model describing the scattering off the residual magnetic 
moment. 

Model and definition of correlation functions — The 
Hamiltonian of the spin-1 single-channel Kondo model 
reads 

H = js imp m+J2 £ ^ c l c ^ (!) 

key 

where S lmp is the impurity spin-1 located at x — and 
c\ a (c ka ) creates (annihilates) a conduction electron with 
momentum k and spin a. The impurity spin is antifcr- 
romagnetically (J > 0) and locally coupled to the elec- 
tron spin density, s(x) — ipl l (x)^Y-ipfj(x), with ip(x) — 
^2 k e lkx Cka- We limit ourselves here to the case of one- 
dimensional conduction electrons, d = 1. A perturbative 
treatment of the coupling J encounters logarithmic low- 
energy divergencies that result in a breakdown of per- 
turbation theory at an energy scale ksTn ~ De~ 1 ^ uJ \ 
where D is an energy cutoff and v is the electron density 
of states. At energies below Tk, the impurity spin gets 
partially screened from S[ mp = 1 to S* mp = | by the 
conduction electrons which, as a consequence, acquire a 



phase shift of Sq = tt/2. However, the residual moment, 
S* mp , still couples ferromagnetically to the electron spin 
density. The ferromagnetic spin-i Kondo model, com- 
plemented by an additional phase shift of So = tt/2 for 
the conduction electrons is thus expected to capture the 
properties of the Hamiltonian (fT]) at low energies. 

Although a ferromagnetic Kondo coupling is 
marginally irrelevant in the renormalization group 
(RG) sense, it results in important corrections to the 
Fermi liquid properties. Upon RG improvement, the 
effective ferromagnetic coupling acquires an energy 
dependence vJ c s(e) = -1/ log (fesT^/e) for e < ksT K 
where one identifies T K ~ Tk- This logarithmic energy 
dependence is at the origin of the singular Fermi liquid 
properties of the underscreened spin-1 Kondo model. 
Calculating the phase shift for the conduction electrons 
in the Born approximation in J ff one obtains, for 
example, in a small magnetic field, <?/is-B <C ksTx, the 
spin-dependent correction to So = n/2, 

Sa= 2{ 1 + <J 2\og(k B T K / gf , B B))- (2) 

The focus of the present study are the waves within 
the Fermi sea that surge around the spin-1 impurity. A 
snapshot of these is provided by the equal-time correlator 
of the impurity spin and the conduction electron spin- 
density at position x, 

x{x,T) = (Si mp s(x)). (3) 

We also consider the integrated correlator x(T) = 
J dx\(x,T). Using rotational invariance of the Hamil- 
tonian ([I]) , the latter can be expressed as an expectation 
value of the z-component of the impurity spin, S? mp , and 
total spin, S? ot = Sf mp + / dxs z (x), 

X (T) = 3(5f mp ^ ot ) - 2. (4) 

This relation is particularly useful for a numerical evalu- 
ation of x(T) with the help of the NRG. 

Method — Generally, the calculation of spatial correla- 
tions is a non-trivial task since most of the methods used 
to investigate Kondo models are not able to reproduce 
correlation functions. One of the most powerful numer- 
ical methods, Wilson's NRG was believed to be unable 
to handle spatial correlations for a long time. Only very 
recently it has been extended to have good spatial resolu- 
tion by one of the authors @ . The approach is similar to 
which was used to handle the two impurity problem with 
NRG [3]. The key idea behind that method is to map 
the single channel problem onto an effectively two chan- 
nel case in order to have good spatial resolution not only 
around the position of the impurity but around another 
freely chosen point x as well. For the technical details 
we refer the reader to Ref. 0. In such a scheme the com- 
putation of the equal time correlation function is rather 
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FIG. 1: (Color online) Panel a: The equal time spin-spin correlation function, x( x > T) = (Simp s(x)) as a function of the distance 
x measured from the S = 1 impurity for different values of the Kondo coupling j = v J. Unlike the situation for the completely 
screened Kondo impurity, \ changes sign as a function of distance, and the electron spin density close to x — (n + is 
aligned, rather than being anticorrelated, with the impurity spin. In panel b, we show the amplitude of the 2fc_F oscillating and 
uniform part of x at T = that are universal functions of x/£k, see Eq. J5). While for short distances both functions decay as 
~ they both obtain logarithmic corrections for x > £k- As shown in panel c, the envelope of the oscillating part crosses 
over from ~ x~ l to ~ l/[x\n{x / ^k)] at around the Kondo coherence length, £k- Any finite temperature introduces another 
length scale, £t = Kvf/ '(UbT) at which the envelope function crosses over to an exponential decay. The uniform part decays 
faster beyond the Kondo coherence length: it tends to zero as ~ In 2 (x/^k)]- For better visualization, we plot l/[xx(x)] as 
a function of x/£k in panel d. 



simple: x( x iT) appears to be a static thermodynamic 
quantity which can be evaluated with a high precision. 

Results The numerically evaluated correlator 
x(x,T) for different Kondo couplings J (and one- 
dimensional conduction electrons) is shown in Fig. [T^l. It 
exhibits characteristic 2kp oscillations, where kp is the 
Fermi momentum. Very close to the impurity, xkp <C 1, 
as well as at distances corresponding to integer multiples 
of it /kp the spin of the conduction electrons is antiferro- 
magnetically aligned with the impurity spin. In between, 
x{x,T) changes its sign and shows ferromagnetic corre- 
lations. Upon excluding the immediate vicinity of the 
impurity, i.e., distances smaller than n/kp, we find that 
the correlator collapses to an universal curve of the form 

x{x,T)\ x>n/kp = (5) 
X 2kp =- ) cos (2k F x + 25a) + X un — ) 

V sK J-kJ \C,K i-K ) _ 

with two universal functions X2k F and X un . Note that 
this is quite unlike the situation for the Friedel oscilla- 



tions [7| : in that case the non- universal part of the charge 
density oscillations is much more extended in space. The 
additional factor in Eq. (O accounts for the cor- 

rect dimensional units of \(x,T). (Note that s(x) is a 
one-dimensional spin-density.) Moreover, we anticipated 
already the phase shift 5q = it/2 that the electrons attain 
from screening of half of the impurity spin-f . In order to 
identify the two functions Xik F and X un , we indepen- 
dently determined the Kondo temperature Tk- This was 
done by extracting the single particle phase shifts in the 
presence of a local magnetic field B directly from the 
NRG spectra and fitting them to Eq. @. 

The two universal functions X2k F and X un at zero tem- 
perature, T = 0, are shown in Fig. QJ). Both decay alge- 
braically as I j x supplemented by logarithmic corrections 
for large distances, x > ^k- In the presence of a finite 
temperature, the decay changes from algebraic to expo- 
nential for x > where £t = fivp j '(fcgT) is the thermal 
length, as it is shown in panel Fig. [TJ; for the amplitude 
of the oscillating part. 

It is expected that the long-distance behavior of 
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FIG. 2: (Color online) Temperature correction, \(T)— x(0) to 
the integrated spatial spin correlator Q . At low temperatures 
£t < £ic, it follows Eq. ©. 



X(s, T) at J 1 < Tr: is captured by an effective ferro- 
magnetic spin-i Kondo model. Assuming that at large 
distances x 3> ^! the correlator is dominated by the 
unscreened part of the magnetic moment, S? = h, 



X{x, T) ~ {S* s(x)), for x > 



(6) 



we can apply perturbation theory in the effective ferro- 
magnetic coupling J c ff(e) to calculate the asymptotic be- 
havior. The amplitude of the 2Uf oscillating part starts 
in first order and the amplitude of the uniform part in 
second order in J e g(e). Its logarithmic energy depen- 
dence results in characteristic logarithmic corrections to 
power-law decay Q , 



cos (2kpx + 28q) 
log(x/6f) 



for x > £k 



(J) 



Note that the antiferromagnetic correlations of x ob- 
served at x = ?i7r/fcF, with integer n, is an unambigu- 
ous signature of the additional phase shift <5o = tt/2 due 
to the partial screening of half of the impurity spin-1. 
The numerical determination of the coefficient C is hand- 
icapped by the fact that the NRG data only yields re- 
liable information on the spatial spin correlations up to 
x Kk ~ 10 3 such that the asymptotic regime is barely 
reached. In order to compare with Eq. ([7]) nonetheless, 
we employ a two-parameter fit. In addition to the pa- 
rameter C in 0, we allow the Kondo temperature ex- 
tracted from a fit to the leading-logarithmic expression 
of the phase shift (J2J) to vary by an additional constant 
parameter, = AHvf / {kBT^ lascshlft ). Adjusting the 
parameters A and C for a best fit to the long-distance 
tail, we obtain C = 0.07 ± 0.008. 

Within the numerical accuracy, we find the long- 
distance tail of x( x i T) to be described by an universal 
coefficient C in Eq. independent of the Kondo coupling 
J. On the other hand, a renormalization group treatment 
of the effective ferromagnetic Kondo model predicts the 



coefficient C to depend explicitly on the effective bare 
ferromagnetic coupling The numerical finding thus 
indicates that the low-energy model is characterized by 
a universal dimensionlcss ferromagnetic Kondo coupling 
of the order of one. 

This is borne out by considering the integrated corre- 
lator x{T) in Eq. At small but finite temperature, 
T <C Tk, the correction x{T) — x(0) is determined by 
the long-distance tail of the spin correlator T = 0) 
at zero temperature as, effectively, the temperature only 
cuts off the correlation tail at the distance of the thermal 
length £t- From Eq. (0 we thus obtain the asymptotics 



X (T) - X (0) = 



4C 



log(£r/60 



for (8) 



This allows for an alternative determination of the co- 
efficient C by computing xCO with NRG. The result is 
shown in Fig. O It allows a fit to Eq. © over more than 
seven decades in £t/£k yielding a more accurate value 
for C = 0.0666±0.0002. Whereas C is found to be univer- 
sal, the value x(T = 0) depends on the Kondo coupling 
J (not shown). 

To summarize, we have demonstrated that the pro- 
file of the spatial spin correlations, T), around an 
underscreened spin-1 is markedly different from the one 
of a fully screened spin-i. The correlations are more 
pronounced and decay at T = Oasl/xinrf = 1 
with logarithmic corrections, and they change sign as a 
function of distance with ferromagnetic correlations at 
xUf/tt = n + h with integer n. The long-distance be- 
havior is universal and can be explained in terms of an 
effective ferromagnetic Kondo model. 
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